Ab initio Molecular Dynamics Simulations of the Initial Stages of Solid-electrolyte 
Interphase Formation on Lithium Ion Battery Graphitic Anodes 
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The decomposition of ethylene carbonate (EC) during the initial growth of solid-electrolyte inter- 
phase (SEI) films at the solvent-graphitic anode interface is critical to lithium ion battery operations. 
Ab initio molecular dynamics simulations of explicit liquid EC/graphite interfaces are conducted to 
study these electrochemical reactions. We show that carbon edge terminations are crucial at this 
stage, and that achievable experimental conditions can lead to surprisingly fast EC breakdown 
mechanisms, yielding decomposition products seen in experiments but not previously predicted. 



Improving the fundamental scientific understanding of 
lithium ion batteries^— is critical for electric vehicles and 
efficient use of solar and wind energy. A key limita- 
tion in current batteries is their reliance on passivating 
solid electrolyte interphase (SEI) films on graphitic anode 
surfaces*^— Upon first charging of a pristine battery, the 
large negative potential applied to induce Li + intercala- 
tion into graphite decomposes ethylene carbonate (EC, 
Fig. [T]) molecules in the solvent, yielding a self-limiting, 
30-50 nm thick, passivating SEI layer containing Li2C03, 
lithium ethylene dicarbonate ((CH^COsLi^)^"^- 6 - and 
salt decomposition products. C2H4 and CO gases have 
also been detected 7 -! 8 - and shown to come from EC. 9 Sim- 
ilar reactions occur during power cycling when the SEI 
film cracks and graphite is again exposed to EC. 2 If in- 
stead the solvent is pure propylene carbonate (PC), a 
stable SEI film does not materialize^ and the battery 
fails. Our work shows that novel mechanisms for the 
initial stages of SEI-growth at electrode-electrolyte in- 
terfaces can be simulated within time scales accessible 
to ab initio molecular dynamics (AIMD),-- which have 
successfully modelled liquid-solid interfaces^ AIMD is 
likely also applicable to shed light on cosolvent / additives 
which must decompose more readily than EC to alter 
and improve SEI structure, Li + transport, and passivat- 
ing properties*^ 

EC-decomposition mechanisms under electron-rich 
conditions have been proposed (e.g., Refs. 30) an d in- 
vestigated using gas cluster Density Functional Theory 
calculations with and without dielectric continuum ap- 
proximation of the liquid environment J^r— Thus "EC - " , 
coordinated to Li + or otherwise, has been predicted to 
undergo ethylene carbon (CE)-oxygen (Oi) bond cleav- 
age to form a more stable radical anion (Figs, la-b). The 
potential energy barrier involved is at least 0.33 eV i 12 ' 14 
Carbonyl carbon (Cc)-Oi bond- "breaking" (or elonga- 
tion) in the gas phase EC~-Li + complex yields a lower 
barrier, but metastable products^ 

Unlike these previous work, AIMD simulations 
can include explicit liquid state environments and 
EC/graphite interfaces. Unlike classical force field-based 
simulations ) 16 ' 17 AIMD accounts for covalent bond- 
breaking. We apply the VASP code ) 18 ' 19 the Perdew- 



Burke-Ernzerhof (PBE) functional, 20 T-point Brillouin 
zone sampling, 400 eV planewave energy cutoff, tritium 
masses for all protons to allow Born-Oppenhcimer dy- 
namics time steps of 1 fs, and a 10 -6 eV energy conver- 
gence criterion. The hybrid PBEO functional ) 21 ' 22 more 
accurate in many cases, is used for spot checks. A Nose 
thermostat maintains the temperature at T=450 K to 
avoid EC freezing. A uniform neutralizing background 
charge is imposed on systems with net charges.— Non- 
adiabatic quantum effects^ not generally evoked in EC- 
breakdown reactions ; ~ are neglected and will be ex- 
amined in the future. Some gas phase calculations are 
also conducted using the Gaussian suite of programs and 
a 6-311++G(d,p) basis^ 

We first consider liquid EC to motivate subsequent 
breakdown products at the EC-graphite interface. A 
17 ps trajectory with 32 EC molecules and a Li + , but 
no excess e~, is conducted in a (15.24 A) 3 simulation 
cell, corresponding to the experimental 1.32 g/cc EC 
density. 17 Integrating over the first peak of the pair corre- 
lation function (g(r), Fig. QJ) between Li + and O atoms 
yields a solvation number ./Vec = 4.0, in agreement with 
a classical force field predictio n 17 ' 26 but slightly smaller 
than the Raman value of 4.9i2! 

Next we add an excess e~ to three randomly cho- 
sen configurations along the Li + -32 EC trajectory and 
restart AIMD. (Liquid state chemical reactions require 
multiple initial conditions to account for the disorder.) 
The 9.7 kcal/mol Ce-Oi cleavage barrier previously 
predicted-^ implies 10~ 9 s lifetimes at T=450 K if one 
assumes a vibrational prefactor of 10 -13 s -1 . The g(r) 
from one 20 ps AIMD run is indeed consistent with this 
prediction (Fig. QJ). An excess electron is trapped on 
an EC near the Li + . Its Aec averages to 3.3; no bond 
breaking, even of the Cc-Oi typejil is observed. An 
AIMD snapshot confirms that the excess e~ is trapped 
on a EC~-Li + complex^ at the carbonate end (Fig.[TJi). 

In contrast, two other starting configurations lead to 
a new mechanism: a Ce-Oi bond breaks irreversibly in 
two distinct EC molecules outside the Li + first solva- 
tion shell within a surprisingly short 0.8 ps (Fig. [2^,-b). 
Here the excess e~ may have avoided the Li + vicinity be- 
cause the Li + Nec may need to spontaneously decrease 
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FIG. 1: Black, red, blue, and yellow spheres: C, O, H, and 
Li. Orange shading: regions of spin (i.e., excess e~) density 
p 3 > 8xl0" 3 |e|A" 3 . (a) & (b) Isolated EC". Green and vi- 
olet lines depict two breakdown modes; the latter is observed 
only with explicit treatment of liquid EC (see text), (c) g(r) 
between Li + and EC O2 site for systems with and without ex- 
cess e~ . (d) Snapshot of Li + coordinated to intact EC - and 
two other EC molecules (other EC in the liquid omitted). 



to accommodate an electron (Fig. [TJ;), leading to an ki- 
netic barrier. This mechanism may dominate in high 
excess e~ concentration ([e - ]) regions near the anode if 
most Li + there have already complexed with EC~ . While 
the PBE functional we use might slightly underestimate 
barriers compared to hybrid functionals ; 12 i 14 it still pre- 
dicts a barrier of 0.37 eV for isolated EC - and a free 
energy barrier of 0.33 eV when a dielectric continuum 
solvation is added. (See the electronic supporting infor- 
mation, ESI.) Furthermore, gas phase EC" is not ob- 
served to crack within a 7 ps AIMD trajectory. Thus 
explicit treatment of solvent EC molecules, which have 
large dipole moments, appears responsible for the fast 
reaction rate. Given the use of the tritium masses and 
the slightly lower PBE functional reaction barriers, we 
should focus on relative (not absolute) time scales of dif- 
ferent mechanisms, not absolute values. 

The disparity between gas and solution phase rates 
likely arises from the repulsive electron affinity of gas 
phase ECiiS While electronic structure calculations can 
impose a metastable isolated "EC - ," the excess e - lies 
outside the molecule (Fig. [TJd); the Ce-Oi anti-bonding 
orbital is not occupied and bond-breaking is not fa- 
cilitated. In contrast, AIMD/PBE simulations predict 
that, in liquid EC, the excess e~ is initially delocalized 
within one or more EC molecules, unlike in water where 
excess e~ occupies intermolecular spaces stabilized by 
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FIG. 2: (a) In this snapshot, substantial e~ density resides on 
the EC with a Cb-Oi bond about to break, (b) Ce-Oi bond 
breaks in the EC - soon afterwards, (c) & (d) EC~ + e~ 
-> C 2 H 4 + CO^" or 0(C 2 H 4 )OC0 2 -, after 50 fs and 1.0 ps 
respectively. Only a few select EC are shown; the colors are 
as in Fig. [T] except that C and H are grey and white in panels 
(b)-(d). 



hydrogen-bond donors^ 4 - When instantaneous favorable 
molecular geometries localize the e~ on one EC and sub- 
stantially populates orbitals on C and O atoms (Fig. [3^) , 
bond-breaking pathways with rates different from that in 
the gas phase emerge. Using a EC liquid snapshot, we 
confirm that the spatial distribution of excess e~ com- 
puted using PBE and PBE0 functionals are qualitatively 
similar. 25 This behavior depends on the LUMO (low- 
est unoccupied molecular orbital) level; while PBE and 
PBE0 predicts gas phase HOMO (highest occupied MO)- 
LUMO gaps of 6 eV and 8.2 eV respectively, the PBE0 
LUMO energy level is higher by only 0.45 eV. 

Doubling [e _ ] by adding another spin-antiparallel e~ 
to the Fig. [2p liquid configuration yields a C2H4/C0 2- 
pair within 50 fs (Fig. [2b). Starting with the EC~-Li + 
complex in Fig. [TJi and adding an e~ leads to Cc-Oi 
bond cleavage and a 0(C2H 4 )OC0 2_ instead (Fig. [2H; 
see ESI for charge state analyses). This bond was previ- 
ously shown to be the weaker bond when even one elec- 
tron is added to ECf 1 ^ The latter fast 2-e~ mechanism is 
not seen without explicit solvent ^ 

Having re-examined breakdown in liquid EC, we turn 
to interfacial reactions, the main thrust of this work. 
Four layers of Li-intercalated graphite are optimzed in 
a periodically replicated 29.74x14.97x15.06 A 3 cell with 
zig-zag edges exposed (Fig. [3]). The box size preserves 
the density of the liquid region containing 32 EC after 
accounting for the van der Waals radii of the electrode 
atoms. Dangling er-orbitals on edge carbon atoms should 
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FIG. 3: (a) Basal plane LiC 6 plane (CaggLiae); (b)&(c) C=0 
edge (Cig2 48Li72), side and top views; in (b), only the outer- 
most layer of Li chelated to 3 O atoms each are depicted; (d) 
C-OH edge (CmCUsHUsI^). C-H edges (not shown) have 
no Li at edge sites. The stick figures depict graphite C, O, 
and H atoms; color scheme as in Fig. [2] 



FIG. 4: EC breakdown products near the end of AIMD runs. 
Intact EC (not shown) fill the empty spaces, (a) C=0 edge: a 
OC2H4OCO 2 -, and a CO/OC 2 H 4 2 ~ pair (near a Li+ away 
from surface), (b) C=0 (other surface): C03~/C2H4. (c) 
C-OH edge: CO 2- /^^ pair near a Li + away from surface, 
(d) C-OH (other surface): CO/HOC2H4OH, and a H z O from 
disproportionation of 3 C-OH. Color scheme is as in Fig. 



be terminated with hydroxyl (C-OH), quinone (C=0), 
carboxylic acid (COOH) functional groups, and/or pro- 
tons; only the first three are electrochemically active; 2 - 
We perform AIMD simulations on C=0, C-OH, and C- 
H terminated LiCg interfaces, and use static calculations 
on the basal (0001) plane (Fig. [3^), impervious to Li + , 
as a reference. Constant voltage conditions^ are not yet 
feasible with AIMD. Instead, our LiC6 stochiometrySS in 
the slab interior mimics a fully-charged battery anode. 
To some extent our simulation is akin to "immersing" 
SET free LiCg, pre-formed at low voltage, into EC liquid; 
this is analogous to the conditions for SEI growth on Li 
metal, which yield EC-specific SEI similar to those on 
graphitic anodes X 

The optimal geometry when Li occupy all C=0 edge 
sites is depicted in Fig. |3p and |3£. Each Li is chelated 
to two and one O atoms on successive graphite sheets (3 
oxygens total), with edge C=0 bonds tilted out of the 
graphite planes to maximize interaction with Li. The 
mean Li-0 distance is 1.98 A. The Li chemical poten- 
tial (fJ-hi) computed by randomly removing an edge Li 
atom is below -2.4 eV, less than the /iLi=-2.0 eV in the 
basal slab (Fig. [3^), suggesting the number of Li atoms 
is appropriate. For C-OH edge LiCg, when Li occupies 
every sixth edge site chelated with 2 O atoms on suc- 
cessive graphite sheets (4 total, Fig. [3H), ^Li is too high 
at -1.7 eV; too many Li are present. Nevertheless, this 
structure is considered because mixed C=0/C-OH edges 
provide a proton source. The band offsets between these 
model electrodes and isolated EC molecules are more fa- 
vorable for electron transfer than that associated with 



basal plane termination; this is shown in the ESI. 

AIMD simulations confine 32 EC between the C=0 
edge terminated LiCg surfaces (Fig. Et)^ in addition 
to the Li at the C=0 edges, an extra Li + resides in 
the liquid region. Figures Hk-b and Table 1 show that 
e~ transfers from the initially charge-neutral LiC6 into 
EC to form a C 2 H 4 /C0 2_ pair (Ce-Oi cleavage), a 
CO/0(C 2 H 4 )0 2 - pair, and a 0(C 2 H 4 )OC0 2 - (C c -Oi 
cleavages). The last moiety may ultimately crack into 
CO and 0(C 2 H 4 )0 2 ~. The — 2\e\ charge states are con- 
firmed in the ESI. In reactions occuring at electrode sur- 
faces, e~ most likely flows directly to the decomposing 
EC's coordinated to surface Li's without becoming first 
delocalized in the liquid. As such, this mechanism should 
not strongly depend on the EC liquid density or the sim- 
ulation cell size. EC breakdown also occurs in the liq- 
uid region where the dianionic products can bind to a 
Li + . In both cases, the products are consistent with 2- 
e~ addition to liquid EC (Figs.^i & d), except that CO 
formation has not occurred there yet. Future work will 
consider countcrions in the electric double layer, which 
however equilibrate at timescales much slower than the 
observed reactions and are neglected herein. 

For the C-OH edge, both C 2 H 4 + COf - (Fig. Hfc) and 
CO + 0(C 2 H 4 )0 2 " (Fig. SI) products emerge in pi- 
coseconds. The 0(C 2 H 4 )0 2 ~ then extracts H + from the 
electrode to form ethylene glycol, which has been used 
as the chemical precursor (alongside carbonates) to syn- 
thesize (CH 2 C03Li) 2 outside battery settings. 6 To our 
knowledge, CO products, observed in battery charging 
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TABLE I: EC decomposition products in interfacial simula- 
tions, t exclude a 3.5 ps equilibration period during which an 
extra +4 |e| charge is imposed to hinder EC breakdown. In 
Run 2, OC2H4O 2- extracts two protons from C-OH groups 
to form ethylene glycol; a OCHOC2H4OH is formed in Run 3. 



experiments,—^ have not been predicted in calculations 
in the absence of explicit solventsj 

A spin-polarized triplet spin simulation yields qualita- 
tively similar conclusions (Table 1). In contrast, no EC 
decomposition is observed at the H-edge interface in 7 ps. 
The edge dependences are consistent with their known 
electrochemical activities^ but are seldom discussed in 
battery studies because the initial carbon edges become 
masked by SEI formation. 

Our novel predictions suggest the following mechanism 
for initial SEI growth. Fast EC decomposition initiates 
at graphite edge regions rich in oxidized sites £2- Both 
Ce-Oi and Cc-Oi cleavages occur, yielding C0 2 ~ and 
OC2H4O 2 - respectively. OC 2 H 4 2 -, or HOC 2 H 4 OH if 
a proton source exists, potentially reacts with CO2 or 
CO3 - at longer time scales to form (CH2C03Li)2,— a 
main SEI component. Cq-Oi cleavage and OC2H40 2- 
products (Fig. [T^,) are faciliated by Li + and other ionic 
products coordinating to carbon edges. However, slower 



e transfer to the solvent may mandate different decom- 
position mechanisms at later stages of SEI growth^— 

In conclusion, AIMD simulations can yield new in- 
sights concerning ethylene carbonate decomposition in 
electron-rich anode regions. With explicit treatment of 
the EC liquid environment, an excess electron can in- 
duce EC decomposition within AIMD timescales. On 
pristine graphitic electrodes, carbon edge terminations 
strongly affect EC breakdown. C=0 edges give strong 
driving forces and support two EC decomposition path- 
ways, yielding both C2H4 and CO gas and ionic products. 
CO evolution, observed in experiments^ is predicted for 
the first time. Although our reported time scales may 
be somewhat affected by the DFT method used, the re- 
actions pertinent to initial stages of SEI growth are sur- 
prisingly fast, which potentially opens the way for the 
versatile AIMD method to investigate the decomposition 
mechanisms of other solvent or additive molecules^ 
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